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We study a model lipid bilayer composed of a mixture of two incompatible lipid types which 
have a natural tendency to segregate in the absence of membrane fluctuations. The membrane is 
mechanically characterized by a local bending rigidity which varies with the average local lipid 
composition (j>. We show, in the case where k varies weakly with <f), that the effective interaction be- 
tween lipids of the same type can either be everywhere attractive or can have a repulsive component 
at intermediate distances greater than the typical lipid size. When this interaction has a repulsive 
Oh' component, it can prevent macro-phase separation and lead to separation in mesophases with a 

finite domain size. This effect could be relevant to certain experimental and numerical observations 
of mesoscopic domains in such systems. 

: 

PACS numbers: 87.16.-b Subcellular structure and processes 87.16.Dg Membranes, bilayers, and vesicles 
I | 82.70.Uv Surfactants, micellar solutions, vesicles, lamellae, amphiphilic systems 
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V\ , I. INTRODUCTION 

a. 

At the simplest level biological membranes are modeled by homogeneous flexible bilayers of amphiphilic lipid 
I 1 molecules [HQ- However, in many physical and biological situations, these membranes are inhomogeneous on some 
microscopic scale. Indeed four major distinct lipid types are typically present in mammalian cell membranes 
It is natural to ask what may be the role of this homogeneity in the biological context and how it influences the 
mechanical properties of the cell. The interplay between the lipid composition and membrane fluctuations has been 
i i - addressed in many recent studies. The local composition of the membrane will clearly affect its fluctuations and 
local geometry. Indeed, the coupling between membrane fluctuations and local composition is at the origin of the 
budding instability 0,01 seen in certain systems. On the other hand, membrane fluctuations will also influence its 
local composition. In this paper we will examine how the coupling of membrane fluctuations to local composition can 
affect the phase ordering of its component lipids. 

In previous works, the way in which the fluctuation-composition coupling is incorporated into the overall free 
^vq , energy of system falls into two main classes: i) The membrane is composed of a homogeneous lipid background 
t-h ' with added insertions such as trans-membrane proteins and attached polymers, ii) The membrane is modeled as a 
multicomponent system with several lipid types and where the mechanical properties of the system are dependent on 
the relative local concentrations of the various lipid types. 

The insertions considered in models of class i) modify the membrane fluctuations via several different mechanisms. 
First point-like inclusions, such as polymers, exert a pressure distribution on the flexible membrane. This involves 
a coupling of the membrane composition, in this case the density field of the inclusions, to the height h(x) over the 
projected area of the membrane. Another possible coupling is via an imposed boundary condition on the height 
field h at the boundary between the inclusions and the membrane. For example, the contact angle at the boundary 
can be taken to be fixed in order to minimize the hydrophobic free energy of the insertion. This is an example of 
a hard constraint. Alternatively one can introduce a general coupling tensor, related to the orientational degrees of 
freedom of the inclusions, to the local strain tensor V»Vj/i, from which the curvature tensor can be extracted. This 
then corresponds to an energetic term which induces a preferred local curvature. In the literature several types of 
inclusions are considered: circular [E0,H], elliptic more general 0,0,0, El embedded inclusions, as well as 
adsorbed cylinders Besides introducing a tendency for a spontaneous local curvature, which breaks the up-down 
symmetry of the system, inclusions may also modify the energy associated wit terms quadratic in the curvature tensor. 
For example isotropic inclusions may modify the local bending and Gaussian rigidities of the membrane. In the case 
of two inclusions one may then explicitly evaluate their effective interaction. To summarize, the density field of the 
inclusions in all these cases is coupled via: (a) h in the case of insertions exerting a pressure, (b) an effective vectorial 
Vj/i coupling in the case of imposed boundary conditions at the inclusions frontier with the membrane, (c) a two 
tensor coupling to VjVj/i when there is a locally preferred curvature tensor and finally (d) a coupling to ViVj/lVfcVj/i, 
when the local bending and Gaussian rigidities are modified by the inclusions and also when non-isotropic effects are 
present. The above are the most physically relevant couplings up to quadratic order and consequently are the most 
significant in systems where the height fluctuations are relatively small. 

In this case of models of type ii), the variation of the elastic properties of the membrane is more continuous 
than in the case of inclusions. If one neglects the possibility of non-isotropic effects, the most natural parameters 
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which will vary with local lipid composition are the bending rigidity k, the Gaussian rigidity k and the spontaneous 
local curvature c. For instance a concentration dependent spontaneous curvature is considered in [lfil IT^JTR IT^ | . 



Linear perturbations to both the bending rigidity and the spontaneous curvature are studied in [ijj, |2fJ . In [f| , linear 
perturbations to the bending and Gaussian rigidity were considered, the interaction arising in this case is proportional 
to 1 /r 4 and the prefactor is given by the product of coefficients of the linear deviations from the average value of k 
and k. The induced interaction may thus be attractive or repulsive depending on the sign of these coefficients. Of 
course models of type (i), with discrete inclusions, can be described by models of type (ii) when it makes sense to 
take a continuum limit for the inclusions. This limit will be valid for inclusion sizes which are comparable to the 
microscopic length scale of the membrane, that is to say the lateral lipid size. 

All of the studies mentioned above assume zero surface tension. However, the study of membranes under tension 
sheds light upon the physics of biological membranes which are not truly at equilibrium but under external constraints 
or perturbations. The surface tension can be due to electrostatic interactions with the aqueous solvent or due to the 
presence of molecular protrusions. Furthermore, the external action of laser tweezers on a vesicle attracts phospholipids 
an puts the membrane under tension. This leads to interesting phenomena such as pearling instabilities [21|. It has 
been theoretically shown that the presence of surface tension can induced a repulsive interaction between inclusions of 
the same type [HEl. The model used in |22l |23| | is of type i) and is based on a linear coupling of the inclusion to the 
height of the membrane, for example to model the local pressure exerted by an attached polymer. The interaction is 
sensitive to the strength of membrane-inclusion coupling. In this system, the up-down symmetry of the membrane is 
clearly broken by the linear coupling. Indeed, in many biological situations the up-down symmetry of the membrane 
is clearly broken for instance by different compositions in the top and bottom leaves or by the presence of conical 
trans-membrane inclusions. However, it is interesting to ask if the presence of surface tension can also lead to repulsive 
interaction between domains, with similar lipid composition, even when the up-down symmetry is conserved. 

The physics of phase separation may play an important role in biological systems. It has been experimentally 
shown that erythrocyte membranes which contain many different lipid types form immiscible two dimensional liq- 
uids, which are very close to the miscibility critical point [24[. The resulting thermodynamic forces can affect the 
mechanical properties of the membrane and in particular its shape. However, in turn the fluctuations will also affect 
the distribution of the components in the membrane. As an example, a long-range fluctuation mediated repulsion 
between inclusions, combined with a short-range van der Waals attraction, could lead to the formation of mesoscopic 
domains (l6ll2£| of the inclusions. It has also been shown that the presence of a surface tension modifies the effective 
interactions between conical inclusions |26j | , inclusions of the same type are always repelled but oppositely orientated 
inclusions interact attractively at long distances and then repel at shorter distances. This is in contrast to the case 
where there is no tension when all interactions are always repulsive. 

In this paper, we consider a two-component bilayer with the up-down symmetry and, in general, with a non- 
zero surface tension. We show that for certain variations of the bending rigidity and the local surface energy (the 
composition independent component of which can be interpreted as a surface tension) with the local composition 
in lipids, a fluctuation induced lipid-lipid repulsive interaction can appear between domains of similar composition. 
This, together with a short-range van der Waals attraction, can induce the formation of mesophases. In the scheme 
of previous models, our model falls into the class of type ii) above and our cumulant expansion method is similar to 
that used in and 0, |2(]] • In our study we add an non-zero surface energy, as in [22I |23| , but where this local 
surface energy fluctuates with the local lipid composition. 

The paper is organized as follows. In Section II we present our field theoretical model. In section III using a 
cumulant expansion for small height fluctuations we calculate the induced interaction, this rather technical section 
may be skipped by a reader interested only in the physical consequences of the calculation. In section IV the general 
physical properties and asymptotic behavior of this effective interaction are discussed. Section V is devoted to a 
description of the results which are compared to previous studies. In addition we suggest a possible experiment where 
the effects predicted here could possibly be seen. 



We consider a model membrane with two lipid types A and B and where the top and bottom leaves have the 
same lipid composition. In the most frequent case, at least whenever van der Waals interactions are dominant, it is 
energetically favorable for lipids of the same type to be adjacent. In this case, we can write down a typical attractive 
energy per site E = x4>a4>b where 4>A and 4>b are the liquid volume fraction of lipid A and B and \ > is a Flory 
parameter related to electronic polarizabilities of both molecules. We will consider a coarse-grained model for a field 
</> related to the local surface fraction of the two lipid types, i.e. <fi = 4>a — <Pb, which in the absence of surface 
fluctuations exhibits a continuous phase transition at sufficiently low temperatures. The theory is then described by 
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the Ginzburg-Landau Hamiltonian |l6j | 



(1) 



which is written in a covariant form which ensures the independence of the energy from the choice of the two 
dimensional coordinate system denoted by x. The parameter J is positive and related to the Flory parameter x, it is 
a ferromagnetic interaction and energetically favors lipids of the same type being next to each other. The potential 
V((f) fixes the two characteristic values of 4> and the global composition via chemical potential like terms. As the 
potential V appears in H simply integrated over the area of the membrane, it can be interpreted as a composition 
dependent contribution to the surface energy of the membrane. Indeed, the constant part of V which is V(0) can 
be interpreted as a surface tension because it is coupled to the total physical area of the membrane J ^/gd 2 x. The 
term V(0) can thus be used as a Lagrange multiplier to fix the physical membrane area. As mentioned above, V will 
have a <j> dependence as in the usual Landau models for phase separating systems. As in standard Landau theory, 
we will assume that V is a single well at high temperature and a double well at low temperature. This means that 
the system on a plane will exhibit a continuous phase transition. At the mean field level, this transition occurs when 
when the mass associated to this field theory, given by Mq = V"{4>q), vanishes (where 4>o is the homogeneous mean 
field solution). This transition exhibits a divergent correlation length and corresponds to a macro-phase separation 
which occurs at a critical temperature T — T c . 

In the above, the metric of the membrane surface is denoted by gtj and in the Monge gauge it is given by 

9ij = Sij + V,//V,// (2) 

where h is the height of the surface above the projection plane whose area we will denote by A. The term g denotes 
the determinant of Qij and is given by 

g = l + (Vh) 2 . (3) 

Hence the Hamiltonian given by Eq. already implicitly includes a coupling between the local composition, as 
encoded by <f>, and the membrane fluctuations, as encoded by h. The interface energy or line tension, which corresponds 
to the term quadratic in the gradient, is written as to ensure covariance; g lJ is the inverse of gij and is given by 

9 l + (V/i) 2 [ 1 

Here we note that the fact that one should use the covariant form of the line tension is often forgotten in the literature. 
To lowest, i.e. quadratic order in the fluctuations h, one has 



H I [<j>,h] = J d 2 x ^(V0) 2 +V(<P) + J d 2 x 



J -{v<t>Y{VhY - ■ vh) 2 + \(vh) 2 vw 



(5) 



We now take into account the elastic energy of the membrane so the total Hamiltonian of the system is given by 

H[4>,h] = H'fah] +H s [c/),h}. (6) 
The Hamiltonian for surface fluctuations will be taken to be 

H s [^h} = \ j d*xK{<t>)(V 2 tif (7) 

which is the simplest Hclfrich Hamiltonian for surface fluctuations lj and correspond, strictly speaking, to the first 
term in an V 2 /i-expansion of the mean curvature f2j. This Hamiltonian corresponds to a bending energy with local 
bending rigidity which depends on the local composition characterized by <f>. The two dimensional membrane system 
is assumed to have no spontaneous curvature and thus has an up-down symmetry. More generally, one could also 
include a composition dependence on the Gaussian rigidity, the contribution coming from this term would then cease 
to be a topological invariant and should strictly be included. 

The effective partition function in the presence of membrane fluctuations is given by 

Z = J d[<j>]d[h]exp(-/3H) (8) 

where = ksT is the thermal energy scale. We recall that A is the projected area of the membrane, the physical 
area of the membrane is denoted by A + AA, where AA is often called the excess area. For typical biological 
membranes, AA/A is small, of the order of a few percent, and we will thus legitimately assume, in the rest of the 
paper, that height fluctuations are small compared to the typical length scale of the system. 
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III. CALCULATION OF THE FLUCTUATION INDUCED INTERACTION 



In this section we explicitly calculate the fluctuation induced interaction to second order in the cumulant expansion. 

In the high temperature regime, lipids form a mixed phase characterized by a homogeneous and uniform composition 
4>q, with fluctuations tp about 4>o- I n an ensemble where the average value of (p is fixed write <f> — <p$ + ip where 
<po = A^ 1 J ePx</>(x). Consequently in this case, we have J d 2 x.ip(x) = 0. By assuming that k and V behave 
continuously around 4>q, we expand the total Hamiltonian iJBJ up to 0(ip 2 ) in the fluctuations. This leads to 



H[<f> Q ,4>,h] = AV^o) + H£[cf> ,4>} + H^[p ,h} + AH[p ,^,h], 



H$[h] = i / d 



n(M(v 2 h) 2 + v(<p )(Vhy 



(9) 
(10) 

(11) 



When the term proportional to V(<po) is included in the surface Hamiltonian Hq , as we have chosen to do above, 
V(4>q) can be interpreted as an effective elastic energy. However, because it is constant, V((/)q) can be interpreted as 
an effective surface tension. The part of the Hamiltonian which we will treat perturbatively is 



+ 5 [^(^o)(V^) 2 +«"(0o)(VV]^ 2 



[V'iM^hf + n'iM^h) 2 ]^ 



(12) 



The scheme of the calculation is just slightly different in the case where the value of <po is allowed to fluctuate but 
nothing intrinsically changes. 

We perform a cumulant expansion in the partition function (JHJ as follows 



Z s exp( 



-pAV{fa)) J dty] exp(-f3H') 



1-0{AH) 1 



ft 2 



((AH) 



2\S 



Z^exp(-f3AV(cp )) / dWexpi-pHesifaip}) 



(13) 
(14) 



where = J d[h}exp(—f3H$) and (0) s = (Z^) 1 J d[h]0 exp(— (3H$). The cumulant expansion at this order is 
clearly exact to 0(ip 2 ). The effective interaction at this order is thus given by 







((AH) 



2\S 



(15) 



where the subscript c indicates that it is the connected part of the correlation function. 

Note that only the first term in the cumulant expansion can lead to a quadratic terms in V?/>, however this term 
can be seen to be zero by the following: 



(~(W) 2 (V/i) 2 -(W-V/0 2 )f = -(Vf) 2 ((VI.) 2 ) c s -VrfV^(V 1 W ] ^ 



= -(v^) 2 ((vh) 2 ) s c - v^v.-v % 2((v/0 2 )f = o, 



(16) 



where we have appealed to the isotropy of the system. There is therefore no renormalization of the coupling J. Also 
in the first term of the cumulant expansion, terms linear in ip cancel by definition of ip (as they are integrated against 
a constant by isotropy) and the remaining terms yield (AH)f = i J d 2 xAl 2 ip 2 where the mass Mi is given by 



Ml = \ [V"(M(^hf) s c + K "(fo)((V 2 h) 2 ) s c ] . 
Again to quadratic order in ip, the second order term in the cumulant expansion yields 







((AH) 



2\S 



d 2 xd 2 y^(x)(7(x - y)^(y) 



(17) 



(18) 
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where 



£/(x-y) = 



r /a 



K ' 2 (0o)((V 2 ^(x)) 2 (V 2 fe(y)) 2 ) c s + 2V'(cf, a )K'(M((W 2 h( X )) 2 (Wh(y)f) s c + V' 2 {4>o)({Vh(*)f (Wh(y)f) 

(19) 

The potential C/(x — y) is non-local and characterizes the induced interaction mediated by height fluctuations. The 
various connected correlation functions above are evaluated as 

<(V 2 Mx)) 2 (V 2 My)) 2 )! = -|(V 4 G(x-y)) 2 



2 \S 



<(VMx)) 2 (VMy)) 2 ) c S = -^]T(V,V,G(x-y)) 



((V 2 Mx)) 2 (VMy)) 2 ): = ^(VV 2 G(x-y)) 2 



Here the Green function G is given by 

G = («(0 O )V 4 - vmv 2 )- 1 

where 



V(<h) 



(G — G r , 



m 



1 k(<M 



(20) 



(21) 



(22) 



The intrinsic length £ is usually in the range 10-100 nm for biological membranes. The Green function Go is 
Go = — 7T- In I j- 1 (L is an arbitrary length) and G m is the Yukawa interaction given by 



In two dimensions one has 



-V 2 G ( x)+m 2 G m (x) =5(x). 



G m (x) = -J--ftr (m|x|), 

Z7T 



where -ftTo(a;) is the Bessel function of the second kind of order 0. Using these results we find that 

U(x) =BJ(x)+«(|x|). 



(23) 



(24) 



(25) 



The first term of the r.h.s. is short-ranged with B = — (to 4 /2/3)[k'(</> )/k(</(>o)] 2 [<5(0) — 2m 2 G m (0)] and needs to be 
regularized via an ultra-violet cut-off A = 2ir/a corresponding to a microscopic length scale a which would be of the 
order of the distance between lipid heads. The second term v of the r.h.s. above is the long-range induced interaction 
and is independent of the ultra-violet cut-off. It is given by 



4 r^' 



v(r) 



2(3 



1 2 



K((f>o) 



v a (u) = -^lK%(u) + 2aK£(u) + a 2 



where u = mlxl = mr and 



a = 



V'(M «(0o) 



(26) 
(27) 

(28) 



IV. PROPERTIES OF THE FLUCTUATION INDUCED INTERACTION 



Here we discuss the features of the interaction derived in section III. First, the strength of this interaction is 
polynomial in m and therefore reducing the local surface or elastic energy and increasing the local bending rigidity 
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reduces considerably the interaction energy. This same reduction however also increases the range of the interaction. 
Secondly the interaction strength is set by kgT, which means that this interaction has an entropic origin. 

If a is positive, we see that the interaction between areas having the same lipid type (with the same sign of ip) is 
always attractive. To see this, we note that in the expression for v(u), Eq. <|2ti|l . the coefficients of a°, a and a 2 are 
functions of u which are positive and monotonically decreasing hence yielding an attractive interaction at all distances. 
In a ferromagnetic analogy where the lipid types are characterized by a field cj> having the values concentrated about 
±1, this corresponds to a long-range ferromagnetic interaction which enhances the short-range one already present. 
However when a is negative, v now has a repulsive component, corresponding to the coefficient of a, which could 
prevent macro-phase separation. 

For large u we have 



Ko{u) ~ \/^; ex P(- 



(29) 



and thus at large distances v behaves as 



v(r) 



k B T 
' 4tt 2 



v(<M 



(30) 



This 1/r 4 interaction, which is always attractive, is typically seen between inclusions of the same type in membranes 
without surface tension and is found in many of the studies discussed in the Introduction. However, this 1 /r 4 attraction 
does not have the the same physical origin as in previous studies because, as can be seen by examining the prefactor, 
it is generated solely by the fluctuations of the surface energy. 
For small r we find 



v(r) 



k B T 
1 8ir 2 



2 2 



C(ln 2 r) 



k B T k'^V'^q) 1 
47r 2 k 2 ((/)o) r 2 



C(ln 2 r) 



(31) 



which is again attractive if a > but is repulsive when a < (or k' (<f>o)V (<fro) < 0). In this last case, the overall 
interaction is somewhat frustrated: it is attractive at very short length scales (of the order of the microscopic length 
scale) due to van der Waals interactions (in our model represented by the local ferromagnetic interaction), together 
with a longer range membrane mediated repulsion over intermediate length scales, before becoming attractive at 
longer length scales. One can can suppose that the occurrence of these attractive and repulsive interactions can 
prevent macroscopic phase separation and lead to mesoscopic domains. 



V. MICROPHASE SEPARATION 



The final effective quadratic Hamiltonian can now be written as 

n cS [fa, V] = I I d 2 * [J(ViP) 2 + M 2 ^j 2 ] + \ I d 2 xd 2 y^(xMx - y)V(y), (32) 



2 J Z1 J 2 

where the corrected mass is given by M 2 = Mq + Alf + B. It is important to note that this mass depends on the 
microscopic cutoff a since M 2 includes the expectation values {(y 2 h) 2 ) s and ((Vft.) 2 ) 5 which diverge and must thus 
be regularized and a similar regularization is needed to evaluate B. As already explained in the Introduction, we 
consider systems such that, in the absence of height fluctuations, when Mq — ► the system exhibits a second order 
phase transition with diverging correlation length £q = \fj /Mq, In Fourier space, we find 

Heff [0o, i>] = -^-^ J d 2 q [Jq 2 + M 2 + w(q)] ^(q)^(-q), (33) 

where q — |q| and the Fourier transform and its inverse are defined by 

/(x) - J ^/(q)exp(iq-x), (34) 

/(q) = / d 2 x/(x)expHq-x). (35) 
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■v(0) 



In the Fourier representation the non-local part of the interaction is given by w(q) = v(q) — v(Q) and M 2 
thus gives effective mass for the theory. 

The stability of the homogeneous solution against phase separation is determined by the lipid-lipid correlation 
function in Fourier space (V'(q)V'( c l / )) = (27r) 2 <5(q + lO^W where the structure factor is 



5(q) = 



Jq 2 + Ml + w(q) 



Defining 



.) 



we find three dimensionless parameters in the structure factor 



S(q) 



2ee 



2e/?J[(e/4) 2 + g 2 ] 



W a (q) 



where q — £q, which are a, e(3J and £ 2 = J/M 2 . The Fourier transform of the dimensionless potential is 



W a {q) = 2?r / udu [Jo(qu) - 1] v a (u), 
Jo 



(36) 



(37) 



(38) 



(39) 



where Jo(x) is the Bessel function of the first kind of order 0. A divergence of the structure factor, Eq. I|38|l at q = 0, 
while 5 remains finite for q ^ signals a macro-phase separation. In this case, since we find 5(0) = £ 2 /(f3J), the 
phase separation occurs when the correlation length £ e — > oo and this corresponds to the case where the induced 
interaction does not change the nature of the transition but only changes where it occurs. We show the behavior of 
W a (q) and S(q) for different values of a in Fig^ When a > 0, W a (q) and its first derivative are always positive, and 
the only maximum of the structure factor occurs at q = 0. 



Wa(q) 

0.15 




100 



q 



FIG. 1: Plot of W a (q) for various values of a: 0.05 (dotted line), —0.1 (broken line) and —0.15 (solid line). 



However, when a < we see that for particular values of e/3J and £/£ e , S(q) reaches a maximum at an intermediate 
values of q and this maximum can even diverge at a non-zero wave- vector q*. In this case, the homogeneous solution 
becomes unstable before £ e — > oo which leads to the formation of mesophases (mesoscopic phase separation) with 
a finite characteristic length scale given by £/q*. Note that at large values of q, we have W a {q) ~ alng but this 
short-range component of the induced interaction is dominated by the short-range van der Waals interaction term 
whose strength is controlled by J, and S(q) ultimately decreases as 1/q 2 for large q. The maximum of the structure 
factor diverges for an intermediate wave-vector q* which is implicitly defined by the two following equalities 



W^(q*) + 4ef3Jq* = 
W a (?) + 2e(3J[q* 2 + (Z/£ e ) 2 ] = 



(40) 
(41) 



i.e. when the parabola — 2e(3J[q 2 + (£,/£ e ) 2 ] is tangent to W a (q). Given the number of parameters in our theory 
the evaluation of a complete phase diagram is unfeasible, however the fundamental question we wish to address is 
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whether there is a macro-phase or micro-phase separation. To do this we can examine the structure factor at the point 
where the q = mode becomes unstable, that is to say where the effective mass M e = 0. This is thus equivalent to 
examining temperatures T c which are critical in the true sense. If the modes q ^ are stable at T c then we expect to 
see the macro-phase separation. However if at T c there is already a mode q ^ which is unstable then a micro-phase 
separation must have already occurred at a temperature T > T c . Thus, without having to specify the full theory, we 
can identify when a macro-phase separation is converted to a micro-phase one due to coupling between membrane 
fluctuations and its composition. 

Since we are interested in the behavior of the structure factor when approaching the macro-phase transition (t e — > 
oo); we calculate the onset of the micro-phase separation given by Eqs. 1|40[1 - (|41[1 for £/£ e = i.e. at the critical 
temperature T c . The result is shown in Fig. [21 the gray region corresponds to the region of the phase diagram 
(e/9J, a) where mesophases appears whereas the white region corresponds to macro phase separation. The solid line 
corresponds to the solution of Eqs. (|40l) ~ (|41[l (the dots are the exact solutions). This "phase diagram" is plotted at 
a fixed lipid composition 4>q and fixed temperature T = T c corresponding to the critical point in the (4>, T) space. 
It is important to note that when the temperature deviates from the critical temperature, the correlation length i e 
becomes finite (but very large), and the gray region delimited by the solid line shrinks. However, in the case where, 
for given lipid types, the parameters e(3J and a lie in the gray region, mesophases appear before the macro-phase 
separation at a temperature T > T c . In the extremal situation where we are far from the macro-phase separation the 
region of the phase diagram corresponding to mesophases disappears completely. 

For parameter values belonging to the gray region of the diagram, the structure factor S(q) diverges before q — at 
a finite value q* . When moving along the solid line starting at the origin, the value of q* decreases until we reach the 
point {ej3J — 3.8 10~ 3 , a — —0.557) where q* = 2.095 which is the smallest value and corresponds to a characteristic 
length scale of 5-50 nm for the mesophases. Then q* increases again when |cv| increases. 

2 r £|3J 



0.005 0.01 0.015 




FIG. 2: Theoretical phase diagram of the bi-lipidic fluctuating membrane in the plane (e/3 J, a) at the onset of the macro-phase 
separation, i.e. at the critical temperature T c . The gray region corresponds to the parameter range for which a phase separation 
occurs at a non-zero wave vector q* leading to the formation of mesophases instead of a macrophase separation. 

On the experimental side, the parameters a and e(3J are not easy to determine. The phenomenological parameter 
J coming from the Landau-Ginzburg theory, is somehow related to the van der Waals attractions between lipids and 
varies roughly as 1/T such that [3 J is fixed for given lipid types and independent on temperature. Although we 
do not know exactly the value of this parameter, we can assume (3 J ~ 1. For lipidic vesicles made of a mixing of 
two very different lipids such as DMPC/cholesterol (a long lipid and a short one), the curvature modulus has been 
experimentally measured ko ~ 50/csT (26fcsT with DMPC alone) and increases when the proportion of cholesterol 
increases up to 250 ksT with 50% of cholesterol [2^,|2jJ. Hence we find e(3 J ~ 0.01 which means that the mesophases 
region of the phase diagram can be experimentally reached in such systems. 

Finally, in the mesophase region, the denominator of S(q) given by Eq. i|3tj|) becomes negative and the calculation 
of the preceding sections, based on quadratic fluctuations of ip is no longer valid (the modes q w q* are unstable). 
The nature of the resulting stable mesophase requires further analysis to determine it and this is beyond the scope of 
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the present paper. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have shown that the coupling of membrane composition via a composition dependence on the local 
surface energy and bending rigidity can alter the phase diagram of a membrane composed of a mixture of different 
lipids. Indeed, depending on the physico-chemical properties of the lipids (for instance by modifying the length of 
their hydrophobic tail), the membrane can exhibits a micro-phase separation leading to the formation of so-called 
mesophases at a temperature T > T c , i.e. before an eventual macro-phase separation. In the previous works where 
surface tension was considered the composition-fluctuation coupling was linear and the up-down symmetry of the 
system thus broken. 

The long-range part of the interaction given by Eq. (|26J) behaves as — 1/r 4 for large r. This interaction has the 
same behavior as that found between inclusions in several models where surface tension is not present. For instance 
in tensionless membranes one finds the effective pairwise interaction 0, EJ 



V c (r) =-k B T6Tr 2 (^ 
V r 



(42) 



between circular inclusions with the up-down symmetry, where ro is the radius of the inclusions. The long-range part 
of the interaction found in Eq. H26|l - (|30[) is also proportional to the thermal energy but it is solely due to fluctuations 
in the surface or elastic energy. 

A few works focused on the effect of the surface tension on fluctuation induced interactions 0, |2^, I23L l30| . 
Calculating the potential between two circular inclusions which locally apply a pressure on the membrane, Evans et 
al. found an interaction which is everywhere repulsive [22T | between inclusions of the same type and is given by 



$(r) 



C1C2 

27T7 



K (mr), 



(43) 



where Ci is related to the force distribution of inclusions i acting on the membrane surface and 7 is the surface 
tension. Here again, this interaction is different from Eq. i|26|) in origin but has some similar features: it is present for 
membranes under tension and is repulsive with a typical range of £ ~ 30 nm for biological membranes. Our model is 
very different, it does not assume any pressure distribution acting on the membrane but relies on the behavior of k(4>) 
and V{4>) close to a liquid-liquid immiscibility critical point. This proximity to a liquid-liquid immiscibility critical 
point in a real biological context is supported by beautiful experiments on mono-layers made of lipids extracted from 
erythrocytes |24j . 

In this study we have seen that the induced interaction only has a repulsive when component when a < 0. 
Qualitatively, this means that the signs of k'(</>o) and V'((f>o) are opposite: when a region is locally enriched for 
instance in lipid A, bending rigidity increased (k'(0o) > 0) whereas the effective surface tension decreases (V'(<po) < 0). 
Let us consider for a moment the mean- field theory where one neglects the fluctuations if) about <pQ. Consider an 
incompressible membrane which is constrained to have a constant projected area, for example a membrane supported 
by a frame. Also let the membrane exchange lipid species with the bulk solution around it |3^ |. The mean-field free 
energy as a function of 4>o is given from Eq. H14[) as 



A 



1 

4^/3 



fcdfcln [n(4> )k 4 + F(</> )fc 2 



(44) 



This mean-field free energy must be regularized by the ultra-violet cut-off A. As the membrane is in a solution 
containing a reservoir of lipid species, <po is not fixed but is thermodynamically selected so as to minimize the mean- 
field free energy. In this case, in our previous treatment we should have thus included a term V'(<po)rp in the expression 
for Hq , however this term can be seen to cancel exactly with the first term of the cumulant expansion, which in this 
case is also now no longer zero. The part of the free energy F* which varies with (f>o is given by 



—A— = V{M + 8^3 



A 2 lnU(<A)) + 



A 2 



k(0o) 



1 



(45) 



The calculation carried out in this paper is valid for small surface fluctuations, a way of ensuring that the fluctuations 
are small is by choosing a very stiff membrane. This can be ensured by taking k(</)) large. The equation minimizing 
F*((j) ) can be expressed as 



V'( 



1 



1 



47T/3 J K((j)o)k 2 



47T/3 



kdk 



K{<i> Q )k* + v{<t> Q ) 



0. 



(46) 
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Now physically we must have V(4>o) > 0, as for any effective surface tension (it is necessary to have m real), this 
result implies that the system will naturally be in the region where a < 0. Now it is straightforward to show, see for 
example at the the mean field level used here that ratio of the excess area to the projected area is given as 

AA - 1 ( A kdk 1 1 M0o)A 2 \ 

In terms of the ratio of the excess to projected area, equation 1461 can now be written as 



AA 



A 2 V((/)o)AA 



In the limit where A A/ A is small using Eq. (|28J) we obtain 

A 2 



7T/3 A 



87r/31/(0 o )' 

Now if we write A = lix j a where a is the microscopic length scale we find that 

a = ^ 

2e a 



(48) 



(49) 



(50) 



where e a is the surface energy of a square of the membrane of linear dimension a, i.e. the average surface energy per 
lipid. Hence we find a < which suggests a scenario to observe the formation of mesophases experimentally. For 
instance one can use a membrane composed by a mixture of DPMC/cholesterol and supported by a frame close at a 
temperature close to T c . 

In a more general context, molecular dynamics [ssl ] and Monte Carlo simulations |34j have shown that the bending 
rigidity has a non-monotonic behavior as a function of the short lipid number fraction x s : it first decreases rapidly for 
small x s and then increases slowly, with a minimum around x s ~ 0.6. These studies suggest that for a two-component 
bilayer made of short and long lipids, the gradient of k(0) and V{<f) could have opposite signs but some tuning may 
be required. In this case the effective interaction will have a repulsive component which could induce mesoscopic 
phase separation. 

The issue of mesophase formation has been discussed in several papers. Taniguchi |l6| has shown in a model with a 
linear coupling of the composition <fi to the mean curvature that near to spherical vesicles with off-critical compositions 
exhibit circular domains that closely resemble patterns observed in red blood cell echinocytosis |24j . 

A similar study has been carried out in different geometries [T^ and the same general phenomena are observed. 
Inspired by the problem of pattern formation of quantum dots at the air-water interface, Sear et al. j2|| have studied 
the effects of a short-range attraction (on top of a shorter range hard core) and long-range repulsion in Monte Carlo 
simulations of two dimensional systems of interacting particles. In their simulations both circular domains and stripes 
were observed as is the case in the experiments. 

Finally, by adding an attractive short-range interaction to the potential Eq. I|43|) . Evans et al. have argued that 
mesophase formation |22| could be induced. Hence, it could explain the formation of caveolae buds from cell membranes 
and their striped texture. The mechanism proposed in this paper of course leads to the same phenomenology in the 
case where the effective potential induced by membrane fluctuations has an intermediate range repulsive component. 
However we do not find any repulsion in the situation where a > which implies some conditions on the membrane 
composition which could perhaps be tested experimentally. 

The model presented in this paper can be generalized by considering lipid distributions without the up-down 
symmetry, i.e. with different compositions in the up and bottom leaves. In this case, one would introduce a composition 
dependent spontaneous curvature, c(4>) in the Hamiltonian. If one assumes that the mixed homogeneous phase has 
no spontaneous curvature then one takes c(0o) = and in this case the correction to the long-range interaction is 

v*(r) = -^^{c'(fo)} 2 K (mr) (51) 

and the mass is renormalized (by a repulsive term) . Hence this correction is attractive and could wipe out the above 
repulsive effect. The two-component membrane could also contain trans-membrane proteins. Despite the fact that 
the repulsive interaction between inclusions described by Evans et al. would appear, it is well known that protein 
aggregation also increases the local lipid composition, as observed in erythrocyte membranes where it induces a 
phospholipid enrichment |35l | . The inclusion of protein- like insertions in this two- lipid model could thus produce quite 
rich behavior and is a line worth pursuing. 
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Our study has predicted that it is possible that a membrane whose fluctuations are impeded exhibits a macro-phase 
separation whereas if it is allowed to fluctuate freely this transition becomes a mesophase separation . In a stack of 
membranes the fluctuations are suppressed by Helfrich forces which are of steric origin. Experimentally, therefore, 
one could prepare a stack of bilayers at a lipid composition where the bilayers within the stack exhibit a macro-phase 
separation. However, according to our predictions, a single membrane could possibly exhibit a mesophase separation 
|37| . Another possibility is that one could try and observe the effect predicted here by using charged membranes and 
then varying their rigidity by changing the bulk solution's salt content |38|. 

We emphasize that, in this paper, we have concentrated on an entirely equilibrium mechanism as a possible explana- 
tion for the formation of mesoscopic domains. However in living cells, out of equilibrium effects are of course important. 
Recently the recycling of lipids between the membrane and cell interior has been put forward as a non-equilibrium 
mechanism for the formation of raft-like structures in active systems [3!j |4£j . 
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